Operational Dynamic Modeling Transcending Quantum and Classical Mechanics 



Denys I. Bondar/' ^'[^jRenan Cabrera/ '[^ Robert R. Lonipay,'^'[^Misha Yu. Ivanov,'*'[^and Herschel A. Rabitz^'j^] 

^Department of Chemistry, Princeton University, Princeton, NJ 08544> USA 
^Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 
'Department of Theoretical Physics, Uzhgorod National University, Uzhgorod 88000, Ukraine 
^Imperial College, London SW7 2BW, United Kingdom 
(Dated: April 2, 2013) 

We introduce a general and systematic theoretical framework for Operational Dynamic Modeling 
(ODM) by combining a kinematic description of a model with the evolution of the dynamical 
average values. The kinematics includes the algebra of the observables and their defined averages. 
The evolution of the average values is drawn in the form of Ehrenfest-like theorems. We show that 
ODM is capable of encompassing wide ranging dynamics from classical non-relativistic mechanics 
to quantum field theory. The generality of ODM should provide a basis for formulating novel 
theories. 
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Introduction. One primary goal in science is to con- 
struct models possessing predictive capability. This en- 
deavor is usually achieved by trial and error, with a pro- 
posed model either subsequently revised or completely 
discarded if its predictions do not agree with experimen- 
tal results. Generally such a process is slow, hence au- 
tomatization has been attempted [TJ [2] . 

In this Letter, we develop a universal and systematic 
theoretical framework for Operational Dynamic Model- 
ing (ODM) based on the evolution of dynamical average 
values. As an illustration of ODM's scope, we infer quan- 
tum, classical, and unified quantum-classical mechanics. 
In order to construct a system's dynamical model, we 
first postulate an associated kinematic description con- 
sisting of two independent components: i) the definition 
of the observables' average, and ii) the algebra of the 
observables. ODM applied to observable data, given in 
the form of Ehrenfest-like theorems [see, e.g., Eq. ([T])], 
returns the dynamical model (see Fig. [l] for a graphical 
summary). The system's kinematic description can also 
be deduced from complementary experiments. For ex- 
ample, if the results of a sequential measurement depend 
on the measurements' order, then the algebra of observ- 
ables must be non-commutative [see comments after Eqs. 
([3]) and (13)]. Limited access to experiments capable 
of firmly establishing the kinematics does not preclude 
hypothesizing plausible kinematic descriptions. Some of 
these hypotheses may be rejected within ODM by reveal- 
ing their incompatibility with observable dynamical data 

In the spirit of ODM, starting from the Ehrenfest the- 
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orems [Eq. ([2])], we will obtain the Schrodinger equa- 
tion if the momentum and coordinate operators obey the 
canonical commutation relation, and the classical Liou- 
ville equation if the momentum and coordinate opera- 
tors commute. To establish a link between quantum and 
classical mechanics, we introduce a generalized algebra 
of observables, incorporating both quantum and classical 
kinematics, that ultimately leads to a unified quantum- 
classical mechanics. Most importantly, we will show that 
ODM is applicable to a wide range of physical models 
from non-relativistic classical mechanics to quantum field 
theories, thus making ODM an important tool for formu- 
lating future models. 

Preparing Dynamical Data. In the current work, we 
present the conceptual and theoretical framework of 
ODM putting aside issues of handling noise contaminated 
experimental data. Assume we have multiple copies of ei- 
ther a quantum or classical system (without loss of gener- 
ality we consider single-particle one-dimensional systems 
throughout) . Suppose we can precisely measure different 
copies of the particle's coordinate x and momentum p 
at times {tk}^^i. Upon performing ideal measurements 
of the coordinate or momentum on the n-th copy, we ex- 
perimentally obtain {xn(tk)} and {pn{tk)}, n = 1, . . . , iV, 
requiring a total of 2KN observations. Time interpola- 
tion of these data points returns the functions a;„ (t) and 
Pn{t). We may then calculate the statistical moments 

WW - ^EtiK(t)]' and WW = fEIiKW]' 
for Z = 1, 2, 3, . . . We make the ansatz, resembling a Tay- 
lor series with coefficients ai, bi, Ck,i, di, ei, and fk^i, that 
the first derivative of x(t) 
satisfy 



[x(t)]i and p{t) = [p(t)]i 



dt 



.(t)=^(a,[x(i)]'+6/[pW]')+ E cMmpit)]\ 



dt 
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For non-dissipative quantum and classical systems, these 
relations reduce to 



m-x{t)^p{t), -p{t)^^U'{x){t), (1) 



where -U'{x){t) = Y.idiW)Y- 

Kinematic Description. Generalizing Schwinger's 
motto "quantum mechanics: symbolism of atomic mea- 
surements" [3], we adapt that any physical model is 
a symbolic representation of the experimental evidence 
supporting it. The mathematical symbolism for this pur- 
pose needs to be considered. A formalism specialized to 
describe a specific class of behavior (e.g., classical me- 
chanics expressed in terms of phase space trajectories) 
can be effective, but it may be unsuitable for connecting 
different classes of phenomena (e.g., unifying quantum 
and classical mechanics). In this case a general and versa- 
tile formalism is preferred. Building a formalism around 
Hilbert space is a suitable candidate for this role. Hilbert 
space is well understood, rich in mathematical structure, 
and convenient for practical computations. 

Consider the postulates: i) The states of a system 
are represented by normalized vectors |\1/) of a com- 
plex Hilbert space, and the observables are given by 
self-adjoint operators acting on this space; ii) The ex- 
pectation value of a measurable A at time t is A{t) = 
{'^(t)\A\'^{t))\ iii) The probability that a measurement 
of an observable A at time t yields A is |(A |^'(t)) | , where 
A\A) = A\A)\ iv) The state space of a composite sys- 
tem is the tensor product of the subsystems' state spaces. 
Having accepted these postulates, the rest - state spaces, 
observables, and the equations of motion - can be de- 
duced directly from observable data. Importantly, these 
axioms are just the well-known quantum mechanical pos- 
tulates with the adjective "quantum" removed, as \'^) is 
a general state encompassing classical and quantum be- 
havior. We will demonstrate below that these postulates 
are sufficient to capture all the features of both quantum 
and classical mechanics as well as the associated hybrid 
mechanics. Equation ([T]) rewritten in terms of the axioms 
becomes 



m-{^it)\x\^{t)) = {^it)\p\^it)), 



Inference of Classical Dynamics. Let x and p be self- 
adjoint operators representing the coordinate and mo- 
mentum observables. The commutation relationship 



(3) 



encapsulates two basic experimental facts of classical 
kinematics: i) the position and momentum can be mea- 
sured simultaneously with arbitrary accuracy, ii) ob- 
served values do not depend on the order of performing 
the measurements. In terms of our axioms, the dynam- 
ical observations of the classical particle's position and 
momentum are summarized in Eq. ([2]). 

We now derive the equation of motion for a classical 
state. The application of the chain rule to Eq. ^ gives 

{d^/dt\ X 1^-) + (*| X \d^/dt) ^ (*| p/m |*) , 
{d^/dt\p\^) + {^p\d^/dt) = (*| - U'{x) \^) , (4) 

into which we substitute a consequence of Stone's theo- 
rem (see Sec. l|l]) 



i\d^{t)ldt) ^L\^{t)) , 



(5) 



and obtain 



im{-^{t)\[L,x]\^(t)) = {^{t)\p\^{t)), 

^{^{t)\[L,p\\'^{t)) = - {^{t)\U'{x)\^{t)) . (6) 

Since Eq. (|6| must be valid for all possible initial states, 
the averaging can be dropped, and we have the system 
of commutator equations for the motion generator L, 



im[L,x] — p, i[L,p\ — —U'{x). 



(7) 



Since p and x commute, the solution L cannot be found 
by simply assuming L = L(x,p) (regarding the defini- 
tion of functions of operators see Sec. |ll|. We add into 
consideration two new operators Xx and Xp such that 



[x, Xx] = [p, Xp] = i, 



(8) 



and Xp van- 



and the other commutators among x, p, X^ 
ish. The need to introduce auxiliary operators arises in 
classical dynamics because all the observables commute; 
hence, the notion of an individual trajectory can be in- 



dt 



(*(i)|j5 !*(<)) = {^{t)\ - U'{x) !*(<)) . (2) troduced (see also Sec. [Vl!I|. Moreover, the choice of the 



Koopman and von Neumann [H [5j pioneered the re- 
casting of classical mechanics in a form similar to quan- 
tum mechanics by introducing classical complex valued 
wave functions and representing associated physical ob- 
servables by means of commuting self-adjoint operators 
(for modern developments and applications see Refs. [BJ- 
\7U\). Our operational formulation is closely related to the 
approach proposed in Ref. [21] and recently successfully 
implemented for quantum state tomography [22l|23j. Re- 
garding developments of other operational approaches see 
Ref. 124 and references therein. 



commutation relationships ( 36 ) is unique. Equation ( 36 ) 
can be considered as an additional axiom. Now we seek 
the generator L in the form L = L(x, Xx,p, Xp). Utiliz- 
ing Theorem [T] from Sec. [ll] we convert the commutator 
equations ([7| into the differential equations 

■'nL'^^{x,Xx,p,Xp) ^P, LxJx,Xx,p,Xp) = -U'{x), (9) 

from which, the generator of classical dynamics L is found 
to be 



L = pXx/m - U'{x)Xp + fix,p), 



(10) 
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where f{x,p) is an arbitrary real- valued function. Equa- 
tions ([5]) , p6l ) , and ( 43 1 represent classical dynamics in 
an abstract form. 

Let us find the equation of motion for \{px\'i>{t)) p 
by rewriting Eq. ([5| in the xp-representation (in which 
X — X, ^ —id/dx, p = p, and \p = —id/dp), 



.d_ 

'dt 



p d 



d 



i-—~iU'{x) — 
m ox op 



f{x,p) 



{px\^{t))^0, 
(11) 



which yields the well known classical Liouville equation 
for the probability distribution in phase-space p{x,p; t) — 
\{px\^{t)}\\ 



_9 



p{x,p\ t) 



m ox op 



p(x,p;t). (12) 



Thus, we have deduced the classical Liouville equation 
along with the Koopman-von Neumann theory from Eq. 
([2]) by assuming that the classical momentum and coor- 
dinate operators commute. 

Inference of Quantum Dynamics. The hallmark of 
quantum kinematics is the canonical commutation rela- 
tion 



[x,p] = ih, 



(13) 



which implies i) the Heisenberg uncertainty principle and 
ii) the order of performing measurements of the coordi- 
nate and momentum does matter [3]. The evolution of 
expectation values of the quantum coordinate and mo- 
mentum is governed by the Ehrenfest theorems ([2]). 

We repeat the algorithm exercised in classical mechan- 
ics above. Substituting the definition of the motion gen- 
erator H obtained from Stone's theorem (see Sec. |l| 

ih\dm(t)/dt) ^ H\^{t)) , (14) 

into Eq. ([2|, we obtain 

im[H,x\ = hp, i[H,p\ = -hU'{x). (15) 

Assuming H — H{x,p) and utilizing Theorem [l] from 
Sec. |lll the commutation relations in Eq. ( [T5| ) reduce 
to mHp{x,p) = p and H'^{x,p) = U'{x). Whence, the 
familiar quantum Hamiltonian readily follows 



H = p^/{2m) + U{x). 



(16) 



Since the Schrodinger equation was derived from the 
Ehrenfest theorem s (|2| ) assuming the canonical com- 
mutation relation ( |13[ ), the presentation suggests that 
the Ehrenfest theorems are more fundamental than the 
Schrodinger equation. 

Unification of Quantum and Classical Mechanics. (For 
a detailed discussion see Sec. Ill see also Fig. [2]) The 
fundamental difference between non-relativistic classical 
and quantum mechanics is that the momentum and co- 
ordinate operators commute in the former case and do 



not commute in the latt er [^51 - l?r] . The operators x, p, 



Xx , and Xp obeying Eq. ( 36 1 form the classical operator 



algebra. The unified quantum-classical operator algebra 
is based on Xq, pq, -d^, and -dp satisfying 



[Xq,Pq] = ihK, [Xq,'dx] = [Pq,'&p] = h 



(17) 



^ K ^ 1, while all the other commutators among Xq, pq, 
"dx, and "dp vanish. The operators i!}x and i3p are simply 
introduced so that the quantum algebra (i.e., k = 1) is 
consistent with the classical algebra. The limit k de- 
fines the quantum-to-classical transition with the quan- 
tum algebra smoothly transforming into the classical one 
as K — )• 0. Since h enters in the time derivative of 



Schodingcr equation ( 14 ) as well as in the commutator 
relationship ( 13 1, the limit /i — > encompasses more than 



the criterion that the coordinate and momentum opera- 
tors must commute in the classical limit. This situation 
motivated the introduction of the parameter k. 

As the first step towards unification of both mechanics, 
we apply ODM to 



m-{^{t)\xq\^it)) = {nt)\p,\nt)), 



dt 



{^it)\pq\^{t)) ^ {^{t)\-U'{Xq)\^it)) , (18) 



and obtain the Hamiltonian 



n = 



2m 



U{Xq) 



F (jjq — hndx , Xq + hndj 



(19) 



such that ih\d^i{t)/dt) = ii |*(t)), where F is an arbi- 
trary real- valued smooth function. Note that no Ehren- 
fest theorems for the observables O = O {xq,Pq) can spec- 
ify the function F because [J^,0] = 0. Hence, the func- 
tion F is experimentally undetectable. We shall utilize 
this freedom by finding an F which enforces that the 
Hamiltonian ( 40 1 smoothly transform to become the Li- 



ouvillian (43) in the classical limit. 



The classical and quantum algebras are isomorphic. 
The quantum operators can be constructed as linear com- 
binations of the classical operators in many ways, e.g., 



Xq = X ~ hKXp/2, Pq = p + hKXx/2, 



X., 



Xr. 



(20) 



In particular, demanding that the quantum operators are 
expressed as linear combinations of the classical ones such 
that 



lim Xq — X, lim pq = p, lim 6x — Xx, 



lim 9r, — Xr,, 



lim H — hL, 



(21) 



identifies the function F as (see Theorems |4] and [5] in Sec. 

ml) 



F{p,x) ^ -p^/{2mK) -U{x)/k + 0{1). (k ^ 0) (22) 
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Keeping the leading term in Eq. (61), we show in Sec 



III that only isomorphism (|45j) is compatible with such 
a function F, which leads to the final expression for the 
unified quantum-classical Hamiltonian, 



2m 



U{x,) 



1 



2mK 



Pq 



1 



—p\x + -U X 
m K 



(23) 



that fulfills conditions (60 1. Theorem [6] in Sec. Ill 



states 

that = for any value of k if and only if ?7 is a 
quadratic polynomial. 

We now demonstrate that the Wigner phase-space rep- 
resentation is a special case of the unified mechanics. 
First rewriting the equation of motion 

ih\d^,{i)ldt)^nq^^^{t)) (24) 

in the xAp-representation (for which x = x^ \x — 
—id/dx, p — id/dXp, and Xp — Xp), then introducing 
new variables u — 



transform Eq. ( 24 ) into 



hKXp/2 and v 



hKXp/2, we 



d ihnf 

ihK- — 

dt 2m 



- U{u) + U{v) 



P. = 0, 




(pec l^'^(t)) ^ j dXpp^ I X 



(25) 



Hence, the wave function (px\^i^{t)) is proportional to 
the celebrated Wigner quasi-probability distribution. 



By only demanding a consistent melding of quantum 
and classical mechanics within ODM, we achieved the 
construction equivalent to the Wigner phase-space for- 
mulation of quantum mechanics. The great attraction of 
the Wigner formalism is due to its smooth and physically 
consistent quantum-to-classical and classical-to-quantum 
transitions |26H33j . Our analysis also points to a unique 
feature of the phase-space formulation: no quantum me- 
chanical representation, but Wigner's, has a "nice" clas- 
sical limit. Moreover, since the Wigner function's dy- 
namical equation is recast in the form of a Schrodinger- 
like equation ( 24 ) , efficient numerical methods for solving 



where Pf^(u,v]t) cx (a; Ap |5'„(i)). Therefore, p^, is the 
density matrix for a quantum system with the Hamilto- 
nian ( 16 ) after substituting h — > hn. Note that k enters 
the equation of motion ( |24[ ) as only a multiplicative con- 
stant renormalizing h. From this perspective, the limit 
K — >■ is indeed equivalent to fi. — > 0. The transition from 
the xXp- to xp-representation results in 



flKXp^^ I ^^^P.^] c»PA„ 



the Schrodinger equation may be applied to propagate 
the Wigner function for conceptual appeal and practical 
utility. 

Future Prospects. ODM was introduced to derive equa- 
tions of motion from the evolution of average values and a 
chosen kinematical description. In Sees. |IVpX[ ODM is 
applied to the canonical quantization rule, the Schwinger 
quantum action principle, the time measuring problem 
in quantum mechanics, quantization in curvilinear coor- 
dinates, as well as classical and quantum field theories. 
Additionally, relativistic classical and quantum mechan- 
ics is also melded within this framework in Ref. |34| . 

Variational principles are at the heart of physics. 
Within their framework, the problem of model genera- 
tion is reduced to finding the correct form of the action 
functional, whose Euler-Lagrange equations govern the 
model's dynamics. However, the action is usually neither 
directly observable nor unique; hence, its construction is 
a subject of debate and can only be justified post factum 
by supplying experimentally verifiable equations of mo- 
tion. More important, there are phenomena beyond the 
scope of variational principles (e.g., dissipation). ODM is 
a theoretical framework free of all these conceptual weak- 
nesses since it operates with observable data recast in the 
form of Ehrenfest-like relations. Hence, the equations of 
motion are no longer axioms but are corollaries of the 
'more fundamental Ehrenfest theorems. 
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with Dmitry Zhdanov are much appreciated. 



Supplemental Material for: "Operational Dynamic Modeling 
Transcending Quantum and Classical Mechanics" 

In the main text of the Letter, we introduced Operational Dynamic Modeling (ODM) as a consistent universal 
theoretical framework for inferring dynamical models from observable data (idealized in this work as noise free). 
To construct a system's model, ODM requires: i) the definition of observables' averaging, ii) the algebra of the 
observables, and iii) observable evolution of the average values (see Fig. [T]). The purpose of this supplemental 
material is to employ this technique to encompass a variety of dynamical models not covered in the main text (see 
the list below). Additionally, we provide a detailed derivation of the unified mechanics in Sec. Ill (see Fig. [2] for the 
roadmap of this derivation). 

Before proceeding further, we wish to clarify a few points. There is a widespread belief that the Ehrenfest theorems 
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FIG. 1: The general scheme of Operational Dynamic Modeling. 



Ehrenfest theorems 



^ 

Stone's theorem 

h 




FIG. 2: The derivation of quantum, classical, and unified mechanics within Operational Dynamic Modeling. 
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cannot shed light on the quantum-to-classical transition. Such claims are partially due to terminology. Ehrenfest [3S] 
derived the following 

m| = {^it)\p\^it)) , I = - U'ix) \^it)) , (26) 

which we will exclusively refer to as "the Ehrenfest theorems". However, the same label is often applied to mean that 
the ccntroid of a narrow wave-packet follows a classical trajectory, i.e., 

™| {^{t)\x \^{t)) = {^it)\p |*(t)) , I {^{t)\p !*(<)) « ((*(t)| X . (27) 



While Eq. (26) is a rigorous mathematical identity [36l|37j, Eq. (27) is an assertion based on a physical approximation 
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I. STONE'S THEOREM 

Stone's theorem [40l]4T] can be stated as follows: If U(t) is a strongly continuous unitary group (e.g., which describes 
the evolution of a system), then there is a unique self-adjoint operator H (the dynamic generator, e.g., Hamiltonian 
or Liouvillian) such that 

^^Uit) I/) = HU{t) I/) , (28) 

for all I/) from the domain of the operator H. The latter equation can also be formally expressed as U{t) — exp{—iHt). 

If the time-independent generator of motion H is a self-adjoint operator, then Stone's theorem guarantees not 
only the existence of unique solutions of the time-independent Schrodinger and Liouville equations, but also the 
conservation of the wave function norms. Regarding the generalization of Stone's theorem to the case of a time- 
dependent Hamiltonian see, e.g., Sec. X.12 of Ref. [12]. 

Physically, Stone's theorem is equivalent to assuming that observables smoothly depend on time. 
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II. NONCOMMUTATIVE ANALYSIS: THE WEYL CALCULUS 



Noncommutative analysis |43H56j is a broad and active field of mathematics with a number of important applications. 
This branch of analysis aims at identifying functions of noncommutative variables and specifying operations with such 
objects. There are many ways of introducing functions of operators; however, the choice of a particular definition is 
a matter of convenience [i^ . 

To make the paper self-consistent, we shall review basic results from the Weyl calculus, which is a popular version 
of noncommuting analysis. Theorem [T] plays a crucial role in the current paper. Even though we prove this result 
within the Weyl calculus, it is valid in more general settings (see, e.g., Ref. [17] and page 63 of Ref. [SD]). 

The starting point is the well known fact that Fourier transforming back and forth docs not change a sufhciently 
smooth function of n-arguments. 



/(Ai 



1 A„) 



1 /■ " 

' •' i=\ 



9=1 



(29) 



Following this observation, we define the function of noncommuting operators within the Weyl calculus as 

1 



1 , • ■ • , ^ri 



(27r)'^ 



exp 



i=\ 



9=1 



/(^Ij • ■ • I in). 



(30) 



where the exponential of an operator is specified by the Taylor expansion. 



exp(i) 



fc=0 



(31) 



The identity 



/t(ii,...,i„) = /(il,...,i1,) 



implies that the function of self-adjoint operators (301 is itself a self-adjoint operator. Moreover, one may demonstrate 
that 

f, (ii,...,i„) lim- [/(ii,...,ifc -I- £,..., i„)-/(ii,...,ifc,...,i„) 



1 

(2^ 



(2 



ir]k exp 

1=1 

1 f " 



L 9=1 



1=1 



d 



exp 



L 9=1 



(2 



1 /■ " 

-y, / l[d^idr,ie^p 

' •' 1=1 



9=1 



/L(^li ■ • ■ :Cn)- 



(32) 



Equation (30 ) defines a one-to-one mapping between a function /('Cii ■ • • j ■?«) ^nd a linear operator f{Ai, . . . , A„). By 
the same token, Eq. ( 32 1 establishes a one-to-one mapping between the derivative of a function and the derivative of 
a linear operator. 

The following theorem is of fundamental importance: 

Theorem 1. Let Ai,...,An be some operators and Ck — [Ak,B], k — l,...,n. If [Ak,Ci] = [B,Ck] =0, k,l = 
1, . . . ,n, then 



[/(ii, . . . , i„), B] = Y.^Ak,B]f^^ (ii, . . . , i„). 



(33) 



k=l 



where f{Ai, . . . ,An) is defined by means of Eq. (30). 



8 



Proof. We introduce A := iEq^i'^glA " ^9) and C [A, B] = iYTq=iVqCq] hence, [1,(7] = [B,C] = 0. From the 
following identity: 



[ii • • • i„, B] = ^ ii • • • Ak-i[Ak,B]Ak+i • • • i„, 



k=l 



we obtain /cCyl'="\ It follows from Eq. (311 that 



d 



[exp(i), B] = Cexp(i) = ^ c,-^ exp(i). 



(34) 



(35) 



□ 



Having substituted this equality into Eq. (30), we finally reach Eq. (33). 

In the context of Maslov calculus [371 SSI HO] , theorem [l] has been extended to a more general case where Ck need 
not commute with A„ and B. In Ref. [57], commutators of the type [f{Ai, . . . , An),g{Bi, . . . , _B„)] were considered. 



III. UNIFICATION OF CLASSICAL AND QUANTUM MECHANICS 

We reiterate that the key difference between classical and quantum mechanics is that the operators of momentum 
and coordinate commute in the former case and do not commute in the latter case |25fl27| . The operators x, p, \x^ 
and Ap obeying the commutation relations 

[x, Aa;] = [p, Ap] = i, [x,p] = [x, Ap] = [j5, A^] = Ap] = 0, (36) 

form the classical algebra of operators. Let us introduce another auxiliary algebra: The operators Xq, Pq, 'd^, and ??p 
form the quantum algebra of operators satisfying 

[Xq,Pq]=ihK, (0 K < 1) [Xg , l}.^] = [pg, dp] = i, (37) 

with all the other commutators vanishing. The operators and 'dp are introduced into the quantum algebra so that 
it resembles the classical algebra. 

The values of k in the domain ^ k ^ 1 defines the quantum-to-classical character because the quantum algebra 
smoothly transforms into the classical one as k — 0. Since h enters in the canonical commutator relationship for 
the quantum coordinate and momentum as well as the time derivative in Schodinger equation, the limit — > 
encompasses more than the criterion that the coordinate and momentum operators must commute in the classical 
limit. This situation motivated the introduction of the additional parameter k. 

To understand the transition from quantum to classical mechanics (see Fig. [2] for the roadmap of the current 
section), we first apply ODM to 

m| {^it)\Xg |*(t)) = {^{t)\Pg |*(t)) , I (^it)\Pg \^ {t)) = {t)\ - U'{Xg) \^ {t)) (38) 

and find the Hamiltonian % = 'H{xg,pg,'dx,'^p) such that ih \d'if{t)/dt) — ii \^{t)). Using theorem[l] we derive the 
system of partial differential equations for the function 

whose general solution reads 



K 



F (^pg - hKdx , Xg + hndj^j , (40) 



where -F is an arbitrary differentiable function of two variables. 

Expression (40) was inferred only from the Ehrenfest theorems (38) for the coordinate and momentum. We seek 
to show that F remains free even if the Ehrenfest theorem is known for another observable O — 0{xg,pg). The 
Ehrenfest theorem for O reads 

thf^ {^it)\ 6 |*(t)) = (*(t)| [6,H] \^{t)) . (41) 
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The contribution from F to this Ehrenfest theorem would be measurable if O did not commute with F. However, 

[Pq,Pq - flK-dx] = [Pq, Xq + tindp] = [Xq,Pq - hndx] = [Xq,Xq + tlKdp] = =^> [O, F] = 0; 



thus, the function F is truly undetectable. 

We first set F to zero and consider the Hamiltonian 



1 



^^Ui^Xq) 



to deduce whether the classical Liouville equation 

I \d'^{t)/dt) = L , L = pK/m - U'{x)Xp + fix,p), 

can be recovered from the Schrodinger equation 

ih\d^{t)ldt) =Hq\^{t)) 



(/ is an arbitrary function) 



(42) 



(43) 



(44) 



as K — > 0. 

The classical and quantum algebras are the same, i.e., they are isomorphic. Indeed, the quantum operators (ig, 
Pq, 'dxi and dp) can be constructed as linear combinations of the classical operators [x, p, Aa;, and Ap) in infinitely 
many ways. Three examples of such realizations are: i) Xq = x, Pq — huXx + (hn — l)Ap, iJa; = -^x + Ap, "dp = p — x; 
ii) Xq = X- hnXp, Pq = p, ^ Ax, ^p = Ap; 



iii) Xq=X-hK\p/2, Pq=p + hKXj./2, l^x = Ax, 'dp = \p 



(45) 



The linear isomorphism between the two algebras stimulates the question: Can Xq and Pq be expressed as linear 



combinations of the classical operators such that lim„_j.o Xq 



lim 



«^oPg = Pi 



negatively answers this question. First, we shall prove a more general statement: 



and limK_j.o Hq = hLl Theorem [2] 



Lemma 1. Assume U'{x) ^ 0. If there exist operators Xq and pq such that they are linear combinations ofx, p, \x, 
Xp, and 



lim Xg — X, 
lim Pq — ap, 
hm Hq = phL, 

K, — yO 



(46) 
(47) 
(48) 



then [xq,Pq] — if3{a + l/a)HK + o{k). 

Proof. The asymptotic symbols o{k) and o(l) are defined with respect to the limit k ^ 0. By the condition of the 
lemma, we set 



dxq 
dx 



dXq ^ dXq 



From Eqs. (46)-(48), 



Bh — — = lim — — 
dXx 

whence, we conclude that 



dXx 



I3h 



r5 + ^""^ A 



Pq 



dpq^ , dpq " ,dpq Opq ~ 

"^x H ^^P + TT^ — A„. 



dx 



dXx 



dp 



dXr, 



P 



lim 



Pq dPq 



m K-^o \ Km dXr k 



dXx 



K-i-O V Km OX j. K OXr 



lim 



dpq _ (3h 
dXx a 



K + 0{k), 



dXq . , 

dxr""^^^- 



(49) 



(50) 



(51) 



By the same token, we derive from Eqs. (46)-(48l that 

f)f f)ff 

= lim ^ -^hU'(x) = lim 
5Ap ''^"aAp 



ap dpq 



\ Km dXp K 



1 8t 



(52) 
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Equations ( 46 1 and ( 47 1 imply 



dxg 
dx 



= 1 + 0(1), 



dXq 

dp 



dp 



a + o(l), 



dPq 

dx 



= 0(1). 



(53) 



Substituting Eqs. (51)-(53l into Eq. (491 and calculating the commutator between Xq and pq by using Eq. (36l, we 
finalize the lemma's prooL □ 

Theorem 2 (The strong version of the no-go theorem). Assume U'{x) ^ 0. There are no operators Xq and Pq such 
that they are linear combinations ofx, p, X^, Xp, and lim^^QXq = x, liniK-^oPg =Pj hm^-i-o = ^L, [xq,pq\ — ihn. 

Proof. If such operators exist, then according to lemma [ij [xq,Pq] = 2ihK + o(k), which contradicts the statement of 
the theorem. □ 

Let us consider the dependence of the wave function on k. Theorem [2] implies the following weaker statement, 
which can also be demonstrated independently: 

Theorem 3 (The weak version of the no-go theorem). Assume U'{x) ^ 0. There are no operators Xq and pq such 
that they are linear combinations ofx, p, Xx, Xp, and 

lim (xq ~ x) |*,(i)) = 0, (54) 

(55) 



lim (pq - p) \^,it)) ^0, 



lim I Hq — hL 



(56) 



and [xq,Pq] = ihK. 

This theorem might seem counterintuitive at first sight. To see that the result is correct, we need to elucidate 
the physical meaning of assumptions (54|-(56). The quasi-classical wave function in the coordinate representation is 
known to be of the form 

$^(j:,t) ^ F{x,t)exp[iS{x,t)/{hK)], (57) 
0. Consider the action of the momentum operator. 



where S{x, t) satisfies the Hamilton-Jacobi equation as n 
p = —ihnd/dx, on this wave function 

p^^ix,t)=^^^^^,ix,t)+OiK) 



lim[p~ ^{x,t)]^,ix,t) = 0, 



(58) 



where ^{x,t) :— dS{x,t)/dx denotes a classical particle's momentum. Equation (58) coincides with condition ( [55| , 
which means that the quantum momentum goes over to the classical momentum in the classical limit. Condition (54) 
implies the same for the coordinate. Nevertheless, one readily demonstrates that 



lim 



KHq - jr{x,t) $«(a;,t) = 0, 



(59) 



where J^{x,t) := ^^(x,t)/(2m) + U{x) is the classical Hamiltonian. Equation (59) contradicts condition (56); thus, 
the statement of theorem [2] is intuitively correct because the quantum Hamiltonian does not approach the Liouvillian 
in the classical limit. 

Now we face the dilemma: Hamiltonian ( 40 ) has been introduced as a generator of motion valid in both the classical 
(k 0) and quantum {n — > 1) cases, and yet the classical limit appears to be inconsistent. Condition (56) merely 
seems to demand that such a generalized generator of motion should become the Liouvillian in the classical limit. But, 
it does not. In fact, condition (56) imposes this restriction on Hamiltonian (42), which is a special case {F = 0) of 
more general Hamiltonian (40). The equality limK^o'H = hL is achievable for certain functions F. Before presenting 
a specific example of this function, let us prove the following two statements: 

Theorem 4. Assume F{p, x) — Q{p) + G{x) and U' {x) ^ 0. // there exist operators Xq, Pq, Ox, and Op such that they 
are linear combinations of x, p, Xx, Xp, and 



lim Xq — X, 



lim Pq 

K — 



lim Ox = X^ 



lim Op = Xp, 



lim H = hL, 



then 



F{p,x) 



V 

2mK 



1 



U{x) + 0{1). (k^O) 



(60) 



(61) 
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Proof. The current proof is similar to the proof of lemma [T] From Eq. (60), we have 



, dL dn 

n — — — am — — 



h— — lim 



Km dXx K dXx ' 



dXx 



, dL dU 
= lim 
dXp '^^0 dXp 

- hU'{x) = lim 



Km dXn 



dXn 



Hk 



dpq 

dX„ 



dOx 
dXx 



dxq d9p 
oXx oXx 



Hk 



dBx^ 

dXr, 



F^{p,x) 



dXq 

dX„ 



Kk 



dX, 



(62) 



(63) 



where F{ and F2 denote the partial derivatives of the function F with respect to the first and second arguments, 
respectively [see Eq. (491 regarding other notations]. Due to the assumption F{p,x) = Q{p) + G{x), the expressions 
under the limits in Eqs. (62) and (63) can be represented as the sum of the term depending on x and the term 
depending on p. For Eqs. ( |62[ ) and (63) to be consistent, the first term must be of o(l) in the case of Eq. (62), and 
the second term must be of o(l) in Eq. (63). Since dOx/dXx = 1 + o(l), dOp/dXp = 1 + o(l), dpq/dXx — o(T), and 
dxq/dXp ~ 0(1), we obtain 



h — = lim 

m K-i-o 



P dpq ^ ^, 
Km dXx 



dpq 
dXx 



Hk 



r(p) = -— +0(1), 
mK 



hU'ix) = lim 



dXq 
dXr, 



Hk 



G\x) ^ — U\x) + 0{1). 



(64) 
(65) 



It can be verified that the derived expressions for Q'{p) and G'{x) indeed simultaneously satisfy Eqs. (62) and (63 1 
Hence, we finally reach Eq. (|6T|). 



□ 

Theorem 5. Assume U'{x) ^ 0. // there exist operators Xq, pq, Ox, and 9p such that they are linear combinations 0/ 
x, p, Xx, Xp, and 



lim x„ 



i, lim Pq = p, lim ( 



A. 



lim 9p = Xp 
At— ^0 



lim — hL, 

K-i-O 



dXq . , 



QPq 

9A„ 



o(k). 



(66) 



the 



F{p,x) = -^--U{x) + 0{l). (k^O) 
2mK K 



(67) 



Proof. Equation (|67j) readily follows from Eqs. (62) and (|63j) after taking into account the following estimates: 
dxq/dXx = o{k), dxq/dXp = 0(1), dOx/dXx = 1 + o(TJ, dOp/dXx = o(l), dpq/dXp = o{k), dpq/dXx = o(l), dOx/dXp = 
0(1), and dOpldXp = 1 + o(l). □ 

Theorems |4] and |5] quite explicitly specify permissible forms of the function F. Thus, the leading order term in Eq. 



( 67 ) shall be taken as the definition of the function F. Consider the Hamitonian 

1 



T~Lqc '■— 



2m 



+Uixq) 



— (pq 

JIK \ 



Hk - 



Pi 



2mK 

1 r 

f - 

K 



hK^^ 



K \ 



q + hK-dp 



U (Xq) — U {Xq + hKOp) 



(68) 



Comparing "Hqc [Eq. (68)] with L [Eq. (43l], we deduce that 



Pq~hK9x/2^p, Ox^Xx, axq + f3hK6p = X, 



(69) 



where a and (3 are unknown constants. Substituting Eqs. (69) into Eq. (68) and requiring 



lim Hqc — hL, 



(70) 
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we conclude that only isomorphism (45) between the classical and quantum algebras is permitted. The quantum 
variables Xq and pq are known as the Bopp operators [58] . Hamiltonian ( 68 1 expressed solely in terms of the classical 
operators reads 



^ - ? 1 A hn ^ 
riqc = —pXx + -U \x~ — Ap 

777 Kj \ A 



(71) 



Hamiltonian (68 1 exactly coincides with Liouvillian (43) for some potentials U. Let us find all such cases: 
Theorem 6. T-Lqc = tiL if and only if the potential U is a quadratic polynomial. 
Proof. The xAp-representation of the classical algebra is 



X — X^ \x 



d 



d 



p = i- 



dx' ^ " Q)^^ 

Equation Hqc — HL written in the xAp-representation leads to 

U{x -a)- U{x + a) = -2aU'{x), a := hKXp/2. 
Fourier transforming this equation with respect to x, we obtain 

[auj — sin (aw)] U{lj) — 0, 



(72) 



(73) 



(74) 



where U{u)) — J dx e '^'^'^U{x)/y/2Tr. Since the equation y = sinj/ has the unique solution y = 0, the non-trivial 
solution of Eq. ( 74 ) must be a distribution with support at the origin, whose most general form reads |59j 



(75) 



n=0 



From the identity (see, e.g., Ref. [60] ) 

a;™5(")(a:) = 



ri\m n\ A("-"')(a;) if TO < n, 

^ ^ (n— m)! ^ ' ' 

{—l)^n\6{x) if m — n, 

if TO > n, 



one derives 



[auj - sin(acj)] S'-^'''+^\oj) = ^ (-l)"a2m+i 



m— 1 
ri-l 



2n+ 1 
2to + 1 



[auj - sin(acj)] S'-^^'^uj) = ^ (_l)™a2m+i 



2n 



2m 



From these equations It follows that the first three terms in expansion ( 75 1 are arbitrary. Moreover, substituting 
expansion ( 75 1 into Eq. ( 74 1 , we obtain 



J2 Cn [au ~ sin (auj)] S^^Hu) = ^ /2„+i (a) 5^^^\lo) + /2„ (a) <5(2«-i)(^) ^ q. 



(76) 



rt=0 



n=0 



where 



m— 1 



p + 2to 
2to + 1 



Equation (76) must be satisfied for all a, then fp{a) = 0, Vp. This condition implies that c„ = 0, ri = 3,4,5, 
because themnctions fp are analytic by construction. Therefore, the most general solution of Eq. ( 74 1 reads 

t/(w) = cq6{uj) + ci5'{uj) + C2S"{uj). 
Finally, we note that the inverse Fourier transform of this function is a quadratic polynomial. 



□ 
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IV. DERIVATION OF THE CANONICAL QUANTIZATION RULE 

The equation of motion for a classical observable / = f{p,x) reads 

d/M=I/,^l, (77) 
where — J^{p,x) is the classical Hamiltonian and |-, •] denotes the Poisson bracket 

We apply ODM to the following equation of motion 

I {^it)\ f |*(t)) = {^{t)\ (f^j \^it)) , (79) 

where |/, J^} denotes a quantum analog of the Poisson bracket, which is to be found. 

From Stone's theorem (see Sec. we conclude that there exits a unique self-adjoint operator H such that 



h \d-^{t)/dt) = H !*(<)). Therefore, we obtain from Eq. (79) 



(*(i)| !*(*)) = ^h{^>(t)\ lf,Jfj |*(t)) . (80) 

The stronger version of this equality gives the celebrated canonical quantization rule 

ir^j^hf,H]. (81) 

in 

V. "STRIPPING" THE SCHWINGER QUANTUM ACTION PRINCIPLE 

The Schwinger quantum action principle [51 1611 - 155] states that a propagator's variation is proportional to the matrix 
element of the quantum action's variation, 

(5(a, t \I3, 0) = (i/h) {a, t\ SS{t) |/3, 0) . (82) 
Consider the case when the action's variation is induced by the system's dynamics, i.e., SS{t) ~ ^^^ijp-St. 

^(a, 1 1/3, 0) ^ I {a, 1 1/3, 0) = {z/h) {a, t\ ^ |/3, 0) . (83) 



According to Stone's theorem (see Sec. |T]), 



Equation ( 83 ) can be rewritten as 



«fi^C>(t) = HU{t), \a, t) = U{t) IV-^) . 



-ihj^ (V-ol uHt) 1/3,0) = u\tf-^ 1/3,0) , 



{rP^\U\t)H\fi,0) = U\t)^ 1/3,0) . (84) 
Since the previous equation is valid for any \ijja) and |/3, 0), we conclude that 

which is an operator analogue of the Hamilton- Jacobi equation. 

Therefore, the Hamilton- Jacobi equation lies behind the Schwinger quantum action principle. Note, however, that 
the Schwinger principle is more general than the Hamilton- Jacobi equation because it holds for any type of variations 
(for further details see Refs. [51 ICTHB5] V 
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VI. MEASURING TIME IN QUANTUM MECHANICS 

Time has a peculiar role in quantum mechanics. First and foremost, it is a key independent variable in dynamical 
equations. However, the definition of the self-adjoint operator representing the time observable is quite challenging, 
and it is still a topic of on-going discussions (see, e.g., reviews [64H66] ). In this section we attempt to address the 
problem of defining the time observable in quantum mechanics from the point of view of ODM. Note that all the 
concepts put forward in this section are also applicable to classical mechanics once the Hamiltonian H is substituted 
by the Liouvillian L. 

Time is physically defined and measured by clocks. We come to know the current time by simply observing the 
position of a clock's pointer. In other words, we obtain the value of time indirectly - by measuring some other 
observable. Assume there is an observable, represented by a self-adjoint operator T, such that its average evolves as 

{•^{t)\f\^{t))^at, (86) 
where a is a real constant. Then, the equality t :— T|^) /ct can be taken as the definition of time through the 



observable T. The form of Eq. ( 86 ) allows for the direct application of ODM to find the operator T, 

|(vl'(i)|t|*(t))=a; 

whence, we obtain the commutator equation for the unknown operator 

i[H,f]^ah. (87) 



The classes of operators H and T obeying the canonical commutation relation (87) are generally quite restrictive. 
Pauli "67^ formally demonstrated that the existence of the self-adjoint time operator canonically conjugate to the 
Hamiltonian implies that both operators possess completely continuous spectra spanning the entire real line. This 
statement is known as the Pauli theorem. However, such a theorem does not withstand a thorough and rigorous 
analysis [SH] and is incorrect in general. Recall that the canonical conjugation of the operators T and H is also a 
theoretical justification for the energy-time uncertainty relation. 

There are other ways to define time measurement in quantum mechanics via ODM. We present the following two 
possibilities: First, assuming that there exist a self-adjoint operator 7i and a real valued function /(•) such that 

(vl/(i)|ti !*(<))= at (*(<)! /(i7)|vl/(t)), (88) 
i.e., at := 71 |^) / (^| f{H) \'^), we readily obtain the commutator equation 

i[H,fi]^ahf{H). (89) 
The second method allows for the observable to be time-dependent. Consider the equation 

{m\f2it)\m) = c^it), (90) 

assuming a(t) is an invertible function, such that the instantaneous value of time can be defined as 

t:=a-^(^{^\f2\^)). (91) 

Then, the equation for the unknown operator 72 (t) reads 

i[H, f2{t)] + hfi{t) = fux'{t). (92) 

VII. QUANTIZATION IN CURVILINEAR COORDINATES 
A. Quantum Mechanics in Curvilinear Coordinates 



In this section we extend ODM to n-dimensional curved spaces as well as to Euclidean spaces in curvilinear 
coordinates, which are important special cases. 
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The classical Hamiltonian of interest is a scalar invariant of the form 



H = ^Pf^ffPi^ + U{X). (93) 



Note that the Einstein summation convention is assumed throughout. The classical Hamilton equations of motion are 

X'^ = -f'P,. (94) 



P -.^--±p9jtlp (95) 



Recall that a non-degenerate coordinate transformation 

X'^^X^(X) (96) 
induces the corresponding momentum transformation 

P-P^=^^-^ (97) 

Direct calculations show that the Poisson brackets calculated with respect to {X, P) obey 

[p^, p,i = Ix^ XI = 0, pj = 6^- (98) 



therefore, the transformations (|96| and (97) are canonical. 

According to ODM, the classical coordinates are replaced by expectation values of the corresponding self-adjoint 
quantum operators. However, extra caution is required to consistently define the operators in curvilinear coordinates. 
The probability density is calculated as 

dV = i:*%l^^dX^dX^ ■ ■ ■ dX" ilfij^gd'^X (99) 

with the essential presence of the weight which can be absorbed as part of the wave function by defining the 
weighted wave function 

^ = .gi/V- (100) 
The scalar expectation value is calculated by means of the weighted integral 



(p) = j rp^^d^x, (101) 

which can be expressed as 

(P) = JxP*g--^F (.g-^) ^d^X = J ^g^^F (g^^) d^X. (102) 

According to this scheme, the scalar weighted operator F and the weighted wave function xjj must transform according 
to the following rules 

i/> i/j = 5^/^, F F = g^/'^Fg-^l^, (103) 
such that the expectation value of the weighted operator reads 

xp*Fipd"X, (104) 

which is an invariant scalar under coordinate transformations. With this in mind, the Ehrenfest theorem applied to 
the first Hamilton equation ( 94 ) can be written as 
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where the second term was added in the right hand side to enforce Hermiticity; whence, 



dt 



d ( 1 



1 „ 1 



while the second Hamilton equation ( 95 1 leads to 

d 



d / 1 



According to Stone's theorem (see Sec. |l]), 



dU 
dx<^ 



d 
at 

the derivatives of the expectation values over time are replaced by commutators, 

^^{x^)='-{[n,x^]), j^{p,)='-{[H,p.]). 

Lifting the averaging leads to the following equations for the unknown quantum Hamiltonian: 



d ( 1 



-9 ^p^g-'g '^g^Pp.g 



^ ' dpa ^2171" 

The latter pair of equations reduces to 
dn d f 1 



-9 ^p-^g'g ' g^Pt^g " , 



[n,p^] 



-Q^A^9 ^P.g^g'g^P.g 



dU 
dx'^' 



d f 1 



dx'^ dx'^ \ 2m 



dpa dpa \2m' 
after postulating the canonical commutation relations 

{x%p^A^8-^K 

Therefore, the Hamiltonian reads 

n = l-g-^/%g^/^g^'^gy%g-y^ + U{x) 



-g "p^g^g g^p^g 



du 

dx^ 



(106) 



(107) 



(108) 



(109) 



(110) 



(111) 



(112) 



(113) 



The problem of constructing consistent quantum Hamiltonians directly in curvilinear coordinates without having to 
rely on the Hamiltonian in Cartesian coordinates was solved by Podolsky [53] ■ In particular he derived Hamiltonian 

( fTi3t . 

Note that the approach presented above based on ODM is equivalent to standard tensor calculus methods. By 
definition, a scalar $ of weight N transforms as 



with the covariant derivative being 



3a; 
dx 



N 



NT"' * 



Similarly, a scalar $ of weight 1/2 transforms as 



such that the covariant derivative reads 



V^* = 1^ - :^r" * = g'/%g-'/^^. 



(114) 



(115) 



(116) 



(117) 
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As a result, the following transformation rule takes place 

V,^^V,^^g'/'—V,^. (118) 
The contraction of a contravariant tensor of weight 1/2 with the covariant derivative is 

VpF'' = ^ + 2^''rV = .g-'/'a^ (.g'/V^J . (119) 
This expression can be used to construct a scalar of weight 1 /2 by contraction, 

V.g^'V,^ = g''^% [g'/'g^'-g'/^d^g''/''^) , (120) 
such that the following transformation rule is obeyed 

V^S^'^V,* ^ V^g^^'V,* = gi/4V^5'"'V,*, (121) 

which ultimately leads us to the conclusion that J ^^g^^^V ^'^cPx is an invariant scalar under coordinate transfor- 
mations. 

In the context of weighted wave functions, the partial derivative is neither invariant under transformations nor a 
scalar. This problem can be partially solved by defining the weighted momentum operator through the covariant 
derivative, 

= -iW^. (122) 

The explicit form of the covariant derivative is problem dependent, e.g., if the wave function is a scalar of weight 1/2, 
then the covariant derivative takes the form 

V^V' = .9'/'5^.9-'/V- (123) 

Even though such a momentum operator is properly defined, it may not necessarily be self-adjoint, and therefore, 
may not possess physical expectation values. Moreover, it is more natural to calculate the expectation values of the 
contravariant components, which are dimensional quantities with fixed physical meanings. 



B. Phase Space Representation in Curvilinear Coordinates 



Let us derive the phase space representation of quantum dynamics in curvilinear coordinates employing ODM. 
Following the recipe in Sec. Ill we begin by extending the quantum algebra with the auxiliary operators 0" and \" 
such that 



[x^',p,]^ihK5^^^, [i^A,]=^5^, [p^,e'']=i5\, [A^,n=0, 

where k is a measure of quantumness/commutativity. 

Note that contrary to Sec. VII A the averaging here does not require the weight 

{F) = [ ^*F^(PXd''P 



(124) 



(125) 



[compare this with Eq. (101)]. According to Stone's theorem (Sec. [T]), the generator of dynamics W in the phase 
space is introduced as 



d 

ih—ip = Wi/'. 
ot 

The generator W must satisfy the following system of equations 

9w law d ( 1 



dW 1 dW 



d f 1 



-Pi^g'^'^PtJ. 



dU 



(126) 

(127) 
(128) 
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The solutions of Eq. ( 128 1 is 



W = 



2mK K V f-y 



(129) 



where / is an arbitrary function. The quantum algebra (124 1 can be realized in terms of the classical operators X^^, 
P^,e^, and A^, 



2 



— A^, 



(130) 



where 



(131) 



The generator W can now be expressed in terms of the classical operators. The function / is specified by requiring 
that the classical limit is recovered as k — >■ 0. Hence, the quantum generator of dynamics in phase space reads 



W = 



1 

2mK, 



2mK 



tlK 

Pn - -r- 



A, 



1 



A,).9'^'^(l + ^e)(p. 



-A,. 



1 



-U X 



-6' 



which can be expanded as 



4m K 

with g^"" = g'"'(X + ^e) and gt" = ^^''(A: - f^O). 
In the A'8-representation 



A„ 



d 



the generator of motion reads 



9T 52 



2mK ae^e"^ 8m^^~ ~ ^+ ' dXt'X'' 8to 

%-' + .9r) 92 ai.g^'^ + .gf) a 1 



Xi^ 



dX-" 



2m 



det^X-" 4m dXt" dQ" k 



u\x-—e]~u\x 



-e 



(132) 
(133) 

(134) 
(135) 



(136) 

(137) 
(138) 



The classical limit is readily calculated from Eq. (1331 



1 

m 



dXt^ 



Applying the XP-representation, 



A„ 



d 



dXt^' 

the classical Liouville equation in curvilinear coordinates 

1 



P^ = P, 



. d 



^g^'•' dp dU dp 



t^^_l„M-pi^ 

dt m^ ^X-'^ 2m''''''dX"dPc,^ dXt'dP^' 



(139) 



(140) 



(141) 



corresponding to Hamiltonian equations ( 94 ) and ( 95 ) , is finally obtained. 
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VIII. CLASSICAL FIELD THEORY 



In the current section, considering the classical Klein-Gordon field, we demonstrate the utility of ODM for classical 
field theories. Conceptually, the case of classical field theories turns out to be similar to the single classical particle 
case. (The Koopman-von Neumann approach has been extended to classical field theories in Refs. [T71 [TOl [20] .) 

The Ehrenfest theorems for the classical Klein-Gordon field read 

i| $ix) = {^it)\7rix) \^{t)) , (142) 

-^{^{t)\7r{x)\^{t)) = {^{t)\^"ix)~f,^{x)Mt)), fi = mc/h, (143) 
c at 

where x can be interpreted as a continuous index and 

mx)J{x')] = [n{x),n{x')] = [${x),n{x')] ^ =^ 

[(^(")(x),(^('")(x')] = [^("Hs;),^^'"^^;')] = [(^^"'(a;),*(")(2;')] =0, Vn,m> 0. (144) 
According to Stone's theorem (see Sec. |l|, we introduce the generator of motion, L, of the state vector, |5'(t)), 

i\d'<i'{t)/dt) ^ cL\^'{t)); (145) 

hence, 

i[L,^ix)]^TT{x), i[L,TTix)] = 4)"{x) - ii^4>{x). (146) 

We shall seek L in the form 

I dx'L(X4x'),\^{x'),${x'),^x'),...J^^\x'),n^"'>{x'),..)j, (147) 

where the auxiliary operators A^(x) and XTr{x) are assumed to obey 

[4>{x), X,f,{x')] ^ iS{x - x'), [7r(x), A^(x')] = i(5(x - a;'), [A0(a;), A^(a;')] = 0. (148) 
Thus, employing theorem [l] from Sec. [llj we get 



L 



dx' [nix')X^ix') + [$"ix') - fi^ix')] A^(x') + (. . . , <^("H^'), ^^"H^O, ••■)}, (149) 



where F = F{. . . , c^*^") (a;), tt^"^ (x), . . .) denotes an arbitrary real functional of derivatives of (j){x) and 7r(a;). 

We shall find the equation of motion for the quantity \{(j){x) tt{x) |5'(t)) P - the probability density for a Klein- 
Gordon field's state being given by (l){x) and 7r{x) at time moment t. We introduce the notation 

^{x) \(f>{x') tt{x')) = <j>{x) |0(x') 7r(x')) , 7t{x) |(/)(x') tt{x')) = 7t{x) \(I>{x') 7t{x')) =^ 

(x) \(j>{x') 7t{x')) ^ {x) \^{x') Tr{x')) , (x) \<P{x') 7r(x')) = ^^"^ (x) |</)(x') 7r(a;')) • (150) 

The functional derivative of F[f{x)] is defined as 

SFlfjx)] _ F[fix)+e6ix-x')]^F[f{x)] 
6 fix') - ~e • 

Whence, 

-f- {F[f{x)]G[f{x)]} = G[f{x)]W^ + Fimf-^]^- (152) 
The operators \(j,{x), Xt:{x), 4>{x), and t:{x) in the 07r-representation read 
A^(.T) = -z^+G(...,0(")(x'),7r("n^'),---), ^-(^) = ~*J^' m = <l^{^), Hx)^7r{x), (153) 
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where G is any real functional. 



n{x) \^{t)) = J dx' |-7r(x')^ - - ^^^{x')] ^ - ^F } {^{x) n{x) . (154) 

Here F has "absorbed" G. 
1 d 



cdt 



This equation is a first order functional partial differential equation. We employ the continuous analogue of the 
method of characteristics to get 



c5t = 



6(t>{x) _ 
it{x) (j)"{x) — ^'^4>(x) 



Sn^x) 



1 S(j){x) 1 SttIx) „ 2,/\ 

7t{x), -—^^^"{x)-fl^(j){x). 



St 



These equations coincide with the classical Klein-Gordon equation 
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x,t) = 0. 



(156) 



(157) 



IX. QUANTUM FIELD THEORY 



Now we shall employ ODM to perform the Bose and Fermi second quantization of the Schrodinger equation. The 
application of ODM to other quantum field theoretic models should be straightforward. 
As before, we start from the Ehrenfest theorems 



I {^{t)\ ^{x) \^{t)) = {^{t)\ - '-Uix)^P{x) + \^{t)) , 

I {^{t)\ i;\x) |*(t)) ^ {^{t)\ '-u{x)i^Hx) - ^^^0^ l*W> • 



(158) 
(159) 



Note that the operator ■0(2;) is not self-adjoint; thus, these Ehrenfest theorems are for complex quantities i^{x) and 
^t(a;). 

Having introduced the generator of motion H by means of Stone's theorem (see Sec. |l| 



ih\d^(t)/dt) ^H\^{t)) , 

with l^'(i)) being an element of a Fock space, we find 

[H,^ix)] = -U{x)i,{x) + [H,i^\x)] = U{x)ijHx) ~ ^ 

Zm ox'^ Zm OX'' 



(160) 



(161) 



A. Bose Quantization 

In the Bose case, we postulate the following commutation relations: 

mx),i:{x')] = [^P^x),^;\x')]^0, [^Jj{x),i:\x')] ^ S{x - x'); 



whence, 



dxjjlx) dtp{x') 
dx ' dx' 



dip^x) drp'fix') 



(162) 



= S'{x' -x), 



dx ' dx 



(163) 
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and seek the generator of dynamics in the form 



dx' ' dx' 



(164) 



Using theorem [T] from Sec. [Tl| and Eqs. (163), we obtain 



^a^t a^ — +'^a*t ^ {^mx) + — — , (165) 



^2„ 



A solution of these equations is 



^U(x)^\x) 



2m dx^ 
2m dx"^ 



^ ^ dj}^ djj ^ cjxy^-^^'^ I C{x)^^\' I 
2to dx dx dx dx 



U{x) + 



dC{x) 
dx 



where C{x) is any real function. Finally, the generator of motion reads 

H = I dx' 



^ dx' 9I/ +Uix')^'i-')^(-') 



where the term that can be represented as the total derivative under the integral was discarded. 

B. Fermi Quantization 

In the Fermi case, the following anti-commutation conditions should be used 

{^(x), ^(x')} = {i^Hx)J\x')} = 0, {4'{x), 4>Hx')} = 5{x - x')- 

whence. 



(166) 



(167) 



(168) 



(169) 



dj'ix') ^ \ _ j d4>{x') d^x)\ _ j d^l^^x') 

dx' '^^''> II dx' ' dx II dx' 



j dti^^x') dtpHxi\ _ 

I dx' ' dx 



= 0, 



j d^Hx') J ' 



dip{x' 



dx 



'^,i,^{a:) \ ^d'{x'-x). 



The crucial difference between the Fermi and Bose cases is the following: 



ipix] 



^I^Hx) 



n 


d4>{x) 


n 


d^\x) 




dx 




dx 



= 0, 



Vn ^ 2, 



(170) 



(171) 



i.e., any power series of ant i- commuting variables terminates after the linear term. Therefore, the generator of motion 
should be sought in the form 



H 



dx' 



dx' dx' 



(172) 



Theorem [T] is not convenient in the Fermi case because it is solely based on commutation relations. However, the 
following identity is a more suitable tool in the case of anti-commuting operators: 



2n 



[il • • • i2„, B] = ^(-l)'=ii • • • Ak-l{Ak, B}Ak+l ■ ■ ■ A2n. 



(173) 



fc=i 



We shall prove this equation by induction. Assuming that Eq. (1731 is correct and utilizing 



[BC,A] = [B,A]C + B[C,A], [BC,A] = -{B, A}C + B{C, A}, 
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we obtain 



[Ai ■ ■ ■ A2nA2n+lA2n+2, B] —[A ■ ■ ■ A2n, B]A2n+lA2n+2 + Ai ■ ■ ■ [^2,1+1 ^2,1+2 , B] 



2n 



: ^(-l)'^il ■■■{Ak,B}--- A2nA2n+lA2n+2 



k = l 



— Ai - ■ ■ A2n{A2n+l,B}A2n+2 + Ai • • • A2nA2n+l{A2n+2-, B} 
2n+2 

= ^(-l)'=il---{ifc,S}...i2„+2. 



fe=l 



Hence, Eq. (|173| is verified 
fow sub 

da2{x) 



Now substituting Eq. (172) into Eqs. (161) and utilizing Eq. (173), we find 



dx 



ai{x) 



- ai{x) 

da*2{x) ' 
dx 



Tpix) -f 

i^\x) 



a2{x) - 03(3;) + 
a2ix) - 02(2;) - 



daajx) 
dx 

da^ix) 



dx 



dx 

djj^jx) 
dx 



as.{x)- 



03(2;) 



dx'^ 

d^i>^{x) 
dx^ 



-Uix)iPix) + 
= C/(a;)^^(a;) 



2m dx^ 
2m dx^ 



Thus, the generator of motion is of the form 



H 



dx' 



di>^x') di>{x') 

2m dx' dx' 



Uix')ip''ix')^lj{x') 



(174) 
(175) 

(176) 



where the term that can be represented as the total derivati ve under the integral was discarded. 
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